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Abstract 

This paper presents an approach to on-line control 
design for aircraft that have suffered either actuator 
failure, missing effector surfaces, surface damage, or 
any combination. The approach is based on a modified 
version of nonlinear dynamic inversion. The approach 
does not require a model of the baseline vehicle 
(effectors at zero deflection), but does require feedback 
of accelerations and effector positions. Implementation 
issues are addressed and the method is demonstrated on 
an advanced tailless aircraft. An experimental 
simulation analysis tool is used to directly evaluate the 
nonlinear system’s stability robustness. 

Introduction 

Reconfiguration in flight concerns the ability of the 
feedback system to redesign itself, in the presence of 
actuator failure, missing surfaces, and wing damage, to 
maintain stability and some acceptable level of handling 
qualities. A reconfigurable control law is generally the 
integration of at least three important elements: 
l)failure detection and isolation to determine which 
surface is no longer useful, 2)on-line parameter 
identification to provide a model of the damaged 
vehicle, and 3)on-line control design which utilizes the 
information from the other two elements to reestablish 
control. This paper deals exclusively with the third 
element: specifically, developing an on-line control 
design methodology that is robust to information errors 
obtained from the other two elements. 

Over the past decade, the literature concerning the 
on-line design portion of reconfigurable controls has 
been dominated by two concepts: dynamic inversion 
(DI) 1 and receding horizon optimal (RHO) 2 control. 
This is not surprising since both can be readily modified 
to handle changing nonlinear dynamics with updates to 
the control’s required onboard model. Also, both 
concepts include provisions for obtaining desired 
closed-loop dynamics (either explicitly or implicitly) to 
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satisfy handling qualities. Unfortunately, both concepts 
tend to be sensitive to onboard model error. 

RHO was used in the Self-Designing Flight 
Control (SDFC) 3 which matured in Lockheed-Martin’s 
approach to the Reconfigurable Control For Tailless 
Fighter Aircraft (RESTORE) program 4 . In SDFC, the 
LQ control stability robustness guarantees of RHO were 
exploited to successfully reconfigure the VISTA/F-16 
on landing subject to a simulated missing left tail 
surface. In RESTORE, the full predictive advantage of 
RHO was utilized to better anticipate the desired plant 
response producing responses that were more 
decoupled and closer to the desired ones. 

Some challenges from SDFC persist for RHO. 
Specifically, the implemented Level 1 models did not 
always result in Level 1 flying qualities. This problem 
was not observed in ground-based simulation. 
Contributing factors could be sensitivity to onboard 
model errors in flight and the lack of guidelines for 
setting or adjusting on-line the algorithm’s cost function 
weightings, penalty terms, and finite time horizon. 
These factors all influence how well the control follows 
the reference model. In contrast, DI avoids the problem 
of manipulating cost functions, but not necessarily the 
problem of sensitivity due to onboard model error. 

The Dl-based, Self-Repairing Flight Control 
System (SRFC) 5 , for example, achieved on-line 
control redesign using onboard models of both the 
nominal and the current/damaged vehicle to detect 
acceleration mismatch. This mismatch drove a pseudo 
inverse allocator to recover the commanded 
acceleration of the nominal flight controller. 

Augmented to an existing F-15 controller, the flight 
experiment demonstrated that this form of dynamic 
inversion was sensitive to model errors not only in the 
current model but in the nominal model which was 
based on a mature aerodynamic database of the F-15. 

A similar issue surfaced when the Intelligent Flight 
Control (IFC) used DI 6,7,8 as an explicit model follower. 
In the IFC structure, dynamic inversion was not applied 
to the actual vehicle but to an onboard model. Plant 
input and desired output reference signals from this 
onboard system were fed to a closed-loop feedback 
system involving the actual vehicle to robustly produce 
the desired response. IFC was not robust enough, 
however, in that it failed to yield the desired 
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longitudinal response on the ACTIVE/F15 due to 
onboard model error produced by the control’s neural 
net. The net was trained on a database for a vehicle 
with square thrust vectoring nozzles. The test vehicle 
had been recently upgraded to round nozzles. 

The issue of model error and dynamic inversion in 
on-line control design was revisited in Boeing’s 
approach to RESTORE 9 . A different explicit model 
following control structure was selected where dynamic 
inversion was applied directly to the vehicle. An on- 
line neural network, based on the work of Kim and 
Calise 10 , adaptively regulated the error in the plant 
inversion to yield the desired rate response of select 
control variables. In truth, the net was a short term 
memory device that modeled the inversion error. If the 
optimal fit (not necessarily obtained on-line) satisfied a 
Lyapunov-based constraint 10 ’ 11,12 , the controller could 
stabilize the vehicle without immediate intervention 
from parameter identification. The optimal fit is 
determined by the net’s structure, i.e. the number/type 
of basis functions and weightings. For the failures 
considered in RESTORE, the number of weightings 
required in the pitch, roll, and yaw channels were 18, 

72, and 72 respectively 13 . For failures not considered, 
no mechanism was offered to adapt the structure, i.e. 
complexity, of the neural net to model different sets of 
inversion errors. This is one approach to desensitize 
dynamic inversion to model uncertainty. Another is the 
approach taken in this paper which reformulates 
dynamic inversion to be less dependent on the onboard 
model. 

The proposed approach to desensitize dynamic 
inversion for on-line control design is an extension of 
Smith’s approach for standard flight control 14 to address 
reconfiguration, various response types, and the 
incorporation of handling qualities. Similar features 
include the absence of an onboard model of the baseline 
aircraft (no control deflections) and the required 
feedback of the accelerations and effector positions. 
Differences include a proposed strategy for dealing with 
nonaffine controls-no inverse mapping via neural nets 
required 10 . Here, a special implementation of 
accelerometers provides the required accelerations with 
considerable redundancy to accommodate sensor 
failure. The proposed control strategy is most closely 
aligned with SRFC but without the required onboard 
models. The proposed control is an implicit model 
following control requiring only the system’s control 
derivatives. Robustness results, using an advanced 
tailless aircraft with multiple innovative control 
effectors 15,16 , will show that the accuracy of these 
control derivatives are not stringent. 

To provide a clear exposition of the method, the 
paper is organized as follows. First, a modified version 


of dynamic inversion is developed for on-line control 
design. Implementation issues associated with 
obtaining the required accelerations, generating the 
vehicles control derivatives, and promoting the desired 
flying qualities follow. Several examples are then 
offered to demonstrate the efficacy of this approach to 
reconfiguration. Lastly, an experimental simulation 
analysis tool is used to directly evaluate the nonlinear 
system’s stability robustness. 

Modified Dynamic Inversion for On-line Control 
Design 

In this section, a general development of the on- 
line control design is offered. Unlike most of the other 
dynamic inversion-based on-line design techniques 
mentioned above, the method pursued here is implicit. 
The desired dynamics of the closed-loop system do not 
reside in some explicit model to be followed but result 
when the feedback loops are closed. Implicit dynamic 
inversion has been used previously in Dl-based flight 
control systems. 1,17,18 Due to the demands placed on the 
on-line control design (i.e., unforeseen rapidly changing 
plant dynamics due to failure occurring possibly during 
complex maneuvering) the form developed here will be 
much less dependent on the control’s onboard model. 

For this discussion, let x denote the state of the 
aircraft whose motion is governed by 

jc = F(x, 8) = f{x) + g(x, 8) (1) 

where the dim(£) > 3. The control objective here is to 
obtain 8 such that three state-dependent control 
variables corresponding to the longitudinal, lateral, and 
directional axes, 

y = bion yiat y<uJ = M*) (2) 

have some desired behavior under both nominal and 
failed conditions. The desired behavior of y to some 

commanded input y c is typically defined by 

y = /*,(* y c ) (3) 

yielding low-order responses that are decoupled along 
axes and compliant with military specifications 19 . 

A dynamic inversion approach to this problem 
involves constructing a control of the form 

8 = yjes) (4) 

to realize some desired control variable rate 
y = y des when substituted into equation 1 subject to 

equation 2. To realize desired closed-loop dynamics of 
equation 3, the control in its implicit form is 

8 = k(x, f des {h(x), y c )) . (5) 

If the input mapping g(x, S) can be expressed as 
g(x, S) = G(x)S , (6) 
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the dynamic inversion control, subject to minimizing 
is 


8 = W(h x G? h x GW(h x G) T \y des - h x f) 


(7) 


where h x = dhjdx and W is a diagonal weighting 
matrix. The dependency on x has been dropped for 
notational convenience. 

The control in equation 7 is clearly dependent on the 
onboard nonlinear model described by / and G . Such 

a model was required for standard flight controls. 17,18 
For on-line control design, the model would have to be 
updated in flight. Due to the model’s form, an update 
requires more information than just the stability and 
control derivatives identified online. 20 

The problem is that the stability derivatives do not 
define the elements of / . Note, the control derivatives 

are the elements of G . The elements of / are defined 


if the on-line parameter identification also identifies a 
bias 3 * , or a model of the form 

jc = Ax + GS + bias (8) 

with 

/ = Ax + bias . (9) 

Here, x is the actual state (not a perturbation) and the 
elements of A and G , are the stability and control 
derivatives, respectively. A key difficulty for dynamic 
inversion is the onboard storage of the nominal bias and 
the required bias corrections which typically depend on 
all the elements of x . 

Another problem with equation 7 for on-line 
control design concerns the system’s actual control 
mapping. The required mapping, GS , is actually the 
nonlinear control mapping g{x,5) rendering the pseudo 
inverse solution in (7) inappropriate. While some have 
proposed neural nets 10 to get around this problem, there 
is a much easier procedure using a modified approach 
to dynamic inversion that eliminates all the problems 
cited. 

To develop the modified form, consider the 
equations of motion as they actually appear with a 
nonlinear control mapping 

F(x,8) = f(x) + g(x,8) = x. (10) 

A standard Taylor series expansion provides the 
following first-order approximation of F(x,S) for x 

and 5 in the neighborhood of [x 0 ,S 0 ] 

F(x, 8) = f(x a ) + g(x a , 8 0 ) + 

JM/w + gU <?)) i x=x 0 ,s=s 0 (*-*„)+ on 

OS-*,,) 


Let x () and 8 0 denote some previous state and control 
from the recent past. The corresponding state rate x Q 


satisfies 

x 0 s f(x 0 ) + 8{x 0 ,S o ). (12) 

Using the standard linear definition, 

A 0 =^(f(x) + g(x,8))\ x=Xo<s=So (13a) 

B °=-^(g( x ^))\x=x 0 ,5=6 0 < 13b > 

with partials evaluated at points (x 0 ,8 0 ) on the 
state/control trajectory, 

F{x,S) = x = x Q + A 0 (x - x Q ) + B 0 AS (14) 


in the neighborhood of [x 0> S 0 \ where AS = 8 - 8 0 . 
With regard to dynamic inversion control, y = /i(x) and 

y = ^x 5 h x (k 0 + \(x - x 0 ) + B 0 AS) . 

djc 

As before, a minimum norm solution provides a 

T -1 

y des - y subject to minimizing AS W AS 

A 8 = W(h x B 0 ) T [h x B 0 W(h x B 0 ) T I”' ( y des - 
(h x x 0 +h x A 0 (x-x 0 ))) 

With a sufficient control update rate, X approaches X Q 
and the control becomes 

A8 = W(h x B 0 ) T [h x B 0 W(h x B 0 ) T ^ \y des - h x x a ) (17) 

with S =S 0 + AS . It is important to note that in this 
formulation x Q is based on the effector position S Q and 
AS is the commanded change to that position in 
response to y des . Hence, sensed positions of the 

effectors are a must. The modified dynamic inversion 
inner loop is shown in figure 1 . Note, the existing S Q 

that corresponds to x Q is taken from the output of the 
actuator. The weighting W consists of the effector’s 
rate limit emphasizing the faster controls in the 
minimum norm solution. 

The desired dynamics defined by equation 3 along 
with equation 17 define the proposed online control 
design method. Note, the control requires the vehicle’s 
control derivatives , B a , either in table lookup form or 

some functional approximation form. The control does 
not require f(x) or g(x,5). Changes in f(x) are 
reflected in the measurement of x a . The online 

identification need only supply the corrected control 
derivatives. The bias and the stability derivatives in 
equation 9 are not required. The problem of applying 
DI to a system with a nonaffine control mapping has 


(15) 
AS so 

(16) 
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Figure 1. Modified DI Inner-Loop 


also been eliminated. Hence, the minimum norm 
solution for the inverse control is viable. Moreover, 
once the failure detection and isolation (FDI) algorithm 
determines which effector has failed, that effector may 
be easily taken off line and the remaining controls 
reallocated by simply setting the corresponding 
weighting to zero. The control does require 
measurements of x Q and S 0 . In the coming section, 

x 0 will be constructed from a special set of 

accelerometer measurements. It is a topic of current 
interest on how to construct clean signals of x Q when 

the accelerometers are subjected to noise, failure, 
biases, and contamination due to structural flexibility. 

Reconfigurable DI Flight Control 

This section includes key equations for assembling 
a reconfigurable control law using the modified form of 
dynamic inversion. Specifically, details are given for 
obtaining x Q and the control derivatives. Two sets of 

control variables will be considered. 

Using a standard x-y-z body axes coordinate frame, 
the equations of motion for a flat-earth, rigid body, 
symmetrical airplane 21 can be expressed as 
ma x ,cg = + % w ~ r b v + 8 sin 0] = X 

ma y cg =m[v + r b u - p b w - g cos 0 sin <f>] = Y 
~ma n cg =m[w+ p b v-q b u- g cos 0 cos <p] = Z 

Pb { xx - r b l xz + % r bV* - V ” Pb% ! xz = L 

% l yy + Pb r b( I xx~ l zz)- r b l xz + Pb 1 xz = M 
ttzz ~ PbU + Pb%Vyy ~ ! xx) + % r b l xz = N < 18 ) 

Here X , Y , and Z represent the total aerodynamic 
and thrust forces (lb), L , M , and N represent the 
total aerodynamic and thrust moments (ft-lb); , lyy , 

, and are the moments and products of inertia 

(slug-ft 2 ); m is the mass (slugs); u , v, and w 
represent the linear velocities (ft/s); u , v , and vv are 
the respective linear accelerations (ft/s 2 ); p b , q b , and 


r b represent the roll, pitch, and yaw rates (rad/s); p b , 
q b , and r b are the respective angular accelerations 
(rad/s 2 ); and a x cg , a y cg . and a ncg are the 

accelerations corresponding to the sum of inertial and 
gravitational forces as sensed by instrumentation 
mounted at the eg. The following auxiliary equations 
must be added to the equations of (18) 

6 - q b oos<p-r b sin </> (19a) 

^ = p b +q b tan 6 sin <f> + r b tan 6 cos <p . (19b) 

The total forces and moments can generally be 
decomposed in the form 

X=X fl +X,+X & (20) 

where X a is the component of axial force due to the 
baseline aerodynamics, X t is the axial force 
component due to engine thrust which includes the 
effects of pitch and yaw thrust vectoring, and X^ is 

the axial force due to aerodynamic effector deflection. 
Let Xft denote that portion of X t due to thrust 
vectoring. 

Consulting equations 18 and 19, let x 0 correspond 
to the state vector 

T 

x = [u. w, q b , 9, V, p b , r b , <p] . (21) 

The linear and angular accelerations will be obtained 
from a special grouping of accelerometers to be 
discussed in the next section. 0 and <j> are obtained 
from equation 19. Some economy can be realized in 
constructing x Q since only those elements multiplying 

the nonzero columns of h x are needed in the control. 
Prior to selecting control variables which determine 
h x , the nonlinear control mapping g(x,S) used to 
define the control derivatives is given. 

Solving for the linear and angular accelerations in 
equation 18, with the decomposition used in equation 
20, and noting the part due exclusively to control yields 
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In equation 22, all angular accelerations are converted 
to units of deg/s 2 . 

In this paper, two sets of control variables are 
considered for reconfiguration. In one set, 

y = h(x) = [q b ,p s ,/3-0.2rJ (23) 


both the forces and moments of the effectors are 
manipulated to produce the desired dynamics. Here, P 


represents the sideslip, and p s and r s represent 

stability axis roll and yaw rates. Force manipulation is 
required due to the ( u , w ) dependence of angle of 
attack, a , in p s and r s and the ( u , v , w ) dependence 

of P . 

In the second set, 

y = h(x) = [q b , p b , r b f (24) 

only the moments are used to produce the desired 
dynamics (not the same desired dynamics as the first 
set). The control resulting from these two sets of 
control variables will be referred to as the ‘Force and 
Moment’ (F&M), and ‘Moment’, (M), approaches to 
reconfiguration, respectively. Details for constructing 
h x can be found in reference 22. 


Accelerometer Measurements 

Determination of the state acceleration vector x 0 

is an important factor in the proposed modified DI 
control law approach. Potential measurement sources 
for the elements of this vector include linear 
accelerometers and numerically differentiated angular 
rates. Both sources have been chosen to get 
independent data. Angular acceleration data is derived 
from a unique implementation of linear accelerometer 
measurements, in addition to the differentiated rates. 
This accelerometer implementation could allow a 
significant amount of sensor redundancy and noise 
reduction, although these aspects have not been 
explored at this time. In addition, linear accelerometers 
are needed for the force and moment (F&M) control 
approach that requires state accelerations. 

In figure 2 three arbitrary sensor groups a t 
(i=l,2,3) are shown, where each sensor group consist of 


n 




-►* 


\ r 3 





a 2 


a 3 


Figure 2. Accelerometer Sensor Group Configuration, 
three orthogonal linear accelerometers ( a xi , , a zi ) 

located a distance r t from the center of gravity (eg). 

The general accelerometer equation, excluding sensor 
noise, bias, and airplane flexibility effects, is 

a i =a cg +mxr i +wxc5'xr i (25) 

where (D is a vector of body-axis angular rates 
( Pb y % , r b )> ® is a vector of body-axis angular 
accelerations (Pb^b^b )> and ^ is a position vector of 
the distances ( , z x ) to the eg. Body-axis angular 
accelerations ( , ^ » ^ ) are obtained by taking the 

difference between any two linear accelerometer sensor 
groups as 

WxAr^j =(ai -a^-axcoxAr^j (26) 

where j represents the distance rj - rj , 

{/, j - 1,2,3; i * One important aspect is that the 

distance between sensor groups, and not the distance to 
the eg, is needed for the angular acceleration solution. 

A total of nine equations can be generated from 
equation 24. In generic form these are 

-r b Ayi_j+q b Azi_j=Ai_j 

V 1 * i-j ~ Pb&i-j = B i-j (2? ) 

“ + Pb^Yi-j = c i-j '■ J ~ 1 - 2 ’ 3 ’ ' * J 

where the right hand side is 

\-j = K., - a x, } r b )Ax ‘-j 

-p b (tb^yi-j-Pb r b^i-j 

Bi-j ={a y ,i -a y j)-Pb‘}b Ax i-j 

, V (28) 

+ (P b +r b )Ayi-j -q b r b Azi_j 

C i—j = ~(. a n,i ~ a n,j ) — P b r b^*i— j 

-Wb^yi-j +(Pb 

with the last equation in (28) written in terms of normal 
acceleration. 

There are multiple solutions for p b , and r b . 

Only five of the nine accelerometers are required with 
the stipulation that at least one accelerometer 
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measurement must come from each of the three sensor 
groups. Sensor redundancy has not been explored in 
this paper, but it appears that many combinations can be 
used for FDI. In addition, by using multiple solutions 
the effects of sensor noise and bias may be reduced by 
some type of averaging technique. 

In this paper, ( ft, , % , r b ) are derived from a least- 

squares solution of 12 equations, 9 from equation (27) 
and 3 equations from the numerically differentiated rate 
gyro signals. All equations are weighted equally, but if 
an FDI system was available equations associated with 
failed accelerometers could be weighted to zero. It is 
assumed that the rate gyros would have a separate FDI 
system. 

The next step is to calculate a cg using equation 25, 
which requires knowledge of . All control designs 
use models that are based upon some fixed eg location, 
so this restriction is not unique. Nine equations are 
solved, 3 equations each for a xcg , a yX g , and a nxg . 

The 3 equations along each axis are averaged to get the 
linear accelerations. 

State accelerations are required for the F&M 
approach and are calculated by rearranging equation 18 
as 


U = a x cg - g sin 9 + vr b - wq b 

(29a) 

v = a y>cg + g sin <j> cos 9 -ur b + wp b 

(29b) 

vv = ~a n cg + g cos <p cos 9 + uq b - vp b 

(29c) 


with units in ft/s 2 . The state acceleration vector i 0 IS 

set up with eight states in the order shown in 
equation 21. 

When the M control approach is used only three 
accelerations are required specifically ( q b% p b , r b ) and 

equation 29 is not used. The linear accelerations at the 
eg are still used, as discussed in the Flying Qualities 
section. 

Control Derivatives 

For this study, complete tables of control 
derivatives have been generated from the aerodynamic 
database. There are at least three reasons for doing this. 

1. The modified dynamic inversion control must 
carry a representation of this control derivative 
database. 

2. Problems with the proposed control redesign 
methodology are separated from potential problems 
with real time system identification methodology. 

3. Tables provide a baseline for the real time system 
identification method. 


The tables are also useful when various uncertainties 
and time delays are injected into the control to quantify 
worst case scenarios regarding the combined FDI and 
parameter identification processes. Here, control 
derivatives are required for both the aerodynamic and 
thrust vectoring controls. Thrust vectoring control 
derivatives result analytically from partials taken of 
functional expressions in reference 22. Only the 
aerodynamic control derivatives will be considered in 
this paper. 

In the typical airplane simulation, force and 
moment aerodynamic data are defined in tables that are 
generated from wind tunnel and flight data. These 
tables, which require interpolation, usually contain 
nondimensional coefficients ( c t ) that are nonlinear 

functions of many airplane variables, for example a , 

(j , Mach, altitude, and S . Representative of the 
aerodynamic forces and moments due to controls, 
consider only the axial force component 

= c&JS (30) 

where q is dynamic pressure (lb/ft 2 ), S is the wing 
reference area (ft 2 ). Aerodynamic force X& has units 
of lb. The term, x , is the nondimensional 

aerodynamic axial force coefficient due to controls. It 
should be noted that c^ x represents a sum from 

several tables including tables of interference 
coefficients representing the effect of one effector 
deflecting upstream from another effector. The c^ x 

terms are also responsible for the nonlinear control 
mapping. Consulting equation 22, the first row of 
dg{x,S)/dS a is constructed from the partials 

dcfa x /dS a . In the simulation, the aerodynamic data 

tables are linearly interpolated. The derivatives, then, 
are constants between the breakpoints of the tables 
(actually up to some € of the breakpoint). The 
generated derivative tables have been constructed with 
this in mind. In testing the on-line control design 
algorithm, the derivative tables are interpolated along 
with the original aerodynamic database tables using the 
simulation’s proprietary interpolation scheme. 

Flying Qualities 

In this section, equation 3 defining the desired 
dynamic behavior for the control variables (CV) is 
considered. Referred to as a ‘command model’, 
equation 3 establishes the flying qualities of the closed- 
loop system. The three control variables and their 
respective command models are considered separately 
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since theoretically the inner DI loop has produced a de- 
coupled integrator block relating y des to y . 

When a slower variable y (d is to be controlled, 

another outer loop can be wrapped around the inner 
loops as illustrated in figure 3. Here the outer loop 
bandwidth is defined as co ol . A nonlinear function, 

NL, is included to relate y 0 i tdes to y c . The reference 

signal is used in the simulation section for comparison 
with the desired system response. 


y<>i,c 


Open-loop 
Command Models 


* 3 ^ 


Reference 

Signals 


(Pol 

— ► 

NL 

y c 

Desired 

y 0 i 

,des 

— ► 

Dynamics 



y 


^ ^des 



Figure 3. Outer loop configuration 


Transfer functions selected for the command 
models are typically found in the military 
specifications 19 . Command models in this paper vary 
with selection of either the force-and moment approach 
or the moment approach. 


In the lateral axis a second outer loop (as illustrated 
in figure 3) is wrapped around the first outer loop for 
p s c to control bank angle <j> during batch simulation 

runs. The two outer loops together create a combined 
transfer function that has the form of equation 33. With 
the second outer-loop frequency defined as co phi , the 

following relationships can be derived as 

tofat = equivalent® 'equivalent (34) 


®phi = 


to equivalent 


2<r 


equivalent 


(35) 


where 0) equimlent and Z equiva ient represent the natural 

frequency and damping ratio of the combined 2 nd -order 
transfer function. These equations are also used in the 
moment approach that follows. Parameter co phi is the 


outer-loop bandwidth for approximating the desired 
bank-angle rate </> des and is defined as 

ides= ®, phifac -*)■ (36) 

A nonlinear equation is inserted into the control loop to 
approximately relate <j> des to p s c as 

Ps,c = fades - (<lb sin <j> + r b cos <p) tan o ) cos <x (37) 

Simulation results will show that these equations give a 
very close approximation to the idealized second-order 
open-loop reference response using co equiva i ent and 


C equivalent * 


Force and Moment (F&M) Command Models 

The longitudinal, lateral, and directional axes 
command-models for the F&M approach are 
respectively 


yion 

^ lontolon( s + to n j on ) 

(31) 



yion,c 

2 2 
s ■*" lon.(Plon s (Plan 


yiat 

_ 

(32) 

yiat,c 

S + to ht 



2 


y dir 

to dir 

(33) 



y dir ,c 

2 2 

$ 2£ dirtodir s to d i r 



Simulations for the F&M approach illustrated in this 
paper include a pitch-rate command ( q c ) for the 
longitudinal axis, a stability-axis roll-rate command 
( p s c ) for the lateral axis, and a linear combination of 

stability-axis yaw-rate ( r s c ) and sideslip-angle ( ) 

commands for the directional axis. The lateral and 
directional axis command models illustrate the mixing 
of linear and angular variables. 


Moment (M) Command Models 

The moment approach only uses angular body- 
rates for the DI loop and only three angular control 
derivatives are required. First-order command models 
are used for p b , q b and r b with transfer functions 

similar to equation 32. A second outer-loop similar to 
the one discussed for the lateral axis above is used for 
all of the control loops. The only difference in the 
lateral loop is the nonlinear equation for body-axis roll 
rate command 

Pb,c - i>des - sin < p +r b cos 0) tan 6 ■ < 38 > 

The longitudinal CV is a and the directional CV is 
f} . Parameters co a and cop are the outer- loop 

bandwidth frequencies for approximating a des and 

P d es respectively as 

<*des =°>a( a c- a ) (39) 

Pdes= C 0 pfa 3 c-P)- (40) 

The nonlinear expressions between a des and q b c and 

between ji des and r bc are 
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gk.cg~ cos<9cos W 

Qb.c des y co$ a CQS p 

p h r tan B 

- tan a tan fijl des + ' 


r - h 
*' c cos a 


*( a >. 


y.cg 


cos a 
costfsin^) 


V cos a cos p 


(41) 


p bc tana. (42) 


In the next section, examples illustrate the 
proposed online control design in various failure 
scenarios for both sets of control variables. Later, an 
experimental simulation tool is used to determine the 
stability robustness of the proposed nonlinear closed- 
loop systems. 


Simulation Results 


Aircraft Model 

A nonlinear simulation using ATLAS* (Aircraft 
Trim, Linearization, And Simulation) that includes a 
highly maneuverable tailless fighter aircraft model with 
innovative control effectors 15,16 has been used for this 
research. Some of the advantages of a tailless aircraft 
are reduced radar signature and reduced weight and 
drag. Analysis of the linearized aircraft model shows 
that it is unstable in the lateral-directional axes at high 
speed. Table 1 shows the longitudinal and lateral 
directional eigenvalues (unstable eigenvalues italicized) 
for three Mach numbers at 25000 ft. 


Table 1. Eigenvalues of Linearized Airplane Model 


Mach 

Longitudinal 

Lateral-Directional 

0.4 

-0.000812 ±j0.0894 

-0.051 13 ±j0. 1609 

-0.5305 ±i 1.667 

-0.2176 ±j0.7351 

0.7 

-0.00328 ± jO.04584 

0.02688 

-0.6621 ± j 1.783 

0.8625 


-0.9073 ± j0.4728 

0.9 

0.1017 

0.009477 

-0.1140 

1.446 

-0.9452 ± j3.556 

-1.303 ±j0.3415 


Eleven control effectors are used in this research, 
but four of these controls are unilateral (operate only in 
one angular direction) leaving nine effective controls. 
These effective controls are: 1) left elevon (LE), 

2) right elevon (RE), 3) symmetric pitch flap (SPF), 

4) all-moving tip (AMT), 5) spoiler-slot-deflector 
(SSD), 6) left outboard leading edge flap (LOLEF), 

7) right outboard leading edge flap (ROLEF), 8) pitch 
vectoring (PV), and 9) yaw vectoring (YV). The AMT 
and the SSD controls are unilateral with control 


* Property of the Flight Controls Branch of Lockheed Martin Tactical 
Aircraft Systems. 


effectors on both the left side (LAMT and LSSD) and 
right side (RAMT and RSSD) of the airplane. 

Second-order actuator dynamics are included in all 
of the simulation runs. These dynamics are composed 
of two real poles with transfer function 


. <»a\ “>a2 

S + CD al s + 0) a2 


(43) 


where the pole values are shown in table 2. 


Table 2. Actuator Dynamics 


Control Effector 

Poles 

(rad/s) 

LE, RE 

40, 100 

SPF 

40, 100 

LAMT, RAMT 

40, 100 

LSSD, RSSD 

40, 100 

LOLEF, ROLEF 

18, 100 

PV 

37,41 

YV 

37,41 


Force and Moment (F&M) Approach 

The first DI control implementation uses a 6-DOF 
force and moment (F&M) approach similar to that used 
in the achievable dynamics research 22 , and discussed in 
a previous section of this paper. Three control variables 
(CV) are defined as shown in equation (23). In this 
vector the longitudinal CV is pitch rate, the lateral CV 
is stability axis roll rate, and the directional CV is a 
linear combination of sideslip and stability axis yaw 
rate. In the lateral axis, a second outer loop is wrapped 
around the first outer loop for p s to control bank angle 
(p during batch simulation runs. 

All command models, corresponding equations and 
variables are discussed in the section on Flying 
Qualities. Values for the command model parameters 
are shown in table 3 with units of rad/s used for all 
frequencies. 


Table 3. Command Model Parameters, F&M Approach 


Parameter 

Value 

Parameter 

Value 

£ Ion 

0.7 

tyat 

2 

<»lon 

5 

£ dir 

0.7 

®n,(m 

3 

“>d<r 

3 

K Ion 

1/3 

® phi 

0.75 


In all of the examples shown below, the controller 
sampling- frequency is 100 hz and the control 
effectiveness matrix is updated from tables every 
iteration. The nonlinear simulation includes the actuator 
dynamics described above, hinge moments, rate limits, 
and position limits. As discussed earlier in this paper, a 
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linear stability matrix A () is not required for this DI 
approach. 

For the baseline example, the airplane was trimmed 
straight and level at 25000 feet and Mach 0.7. A 
directional channel doublet of amplitude 5 degrees was 
commanded during the first 5 seconds. At 7.5 seconds 
into the simulation, a 50 degrees bank angle step was 
commanded, and between 1 1 and 16 seconds a 
5 deg/sec pitch doublet was commanded. In addition, 
at 0.05 second into the simulation, the throttle was 
increased 60° above than the trim value. This large 
increase in throttle along with the downward direction 
of the velocity vector caused Mach number to increase 
from 0.7 to near 0.9 over the 20-second simulation. 

Results are illustrated in the four plots shown in 
figure 4. Solid lines represent the actual control 
variables, long-dashed lines represent input commands, 
and dotted lines represent reference signals from open- 
loop command models. The top plot for the 
longitudinal doublet command shows that the q b 

response almost perfectly duplicates the reference 
signal, except for a little ripple during the p 
movements. The second plot is for the <p response 
which also closely overlays the reference signal. The 
generated error signal is the command p s c that drives 

the p s response shown in the third plot. The slight 
variation between the p s response and its reference 
signal is probably caused by the nonlinear function 
relating $ des and p s c . Finally, the directional doublet 

command, reference signal, and actual response are 
shown in the fourth plot. Tracking results are almost as 
good as the results obtained using the achievable 
dynamics tool 22 that did not include actuator dynamics. 

Figure 5 has the same set of plots for an example 
illustrating a change in aerodynamics resulting from 
two missing surfaces. One second into the simulation 
both the left elevon and the left SSD are declared 
missing, which could happen if part of the left wing is 
damaged. The control system instantaneously knows 
about the missing surfaces and reconfigures by 
excluding these controls. Since the SSD controls are 
unilateral the right SSD is also excluded to prevent the 
control allocator from commanding the missing control. 
Recall that a stability matrix isn’t required since the 
changing forces and moments resulting from the 
revised aerodynamics are completely measured by 
accelerometers. Comparing these results with the 
baseline case in figure 4 shows that tracking is only 
slightly degraded. The system response looks good 
because there is enough control power to compensate 
for the failures. 




30 


<5b 


T3 


0 

-5 




0 5 10 15 20 

Time, sec. 

Figure 4. Baseline force and moment example. 


The ATLAS aerodynamic database incorporates 
changes resulting from missing surfaces. Examination 
of a linear model for this failure condition shows a 
coupled longitudinal and lateral-directional stability 
matrix. Eigenvalues for the double missing surface 
case at Mach 0.7 are shown in table 4. Although the 
longitudinal and lateral-directional eigenvalues cannot 
be isolated, they can be compared with those for the 
Mach 0.7 case in table 1. It appears that the short 
period has changed significantly from -0.6621 ± j 1.783 
in table 1 to -0.8151 ± jl.096 in table 4. The dutch-roll 
(-.9073 ± j0.4728 in table 1) has a smaller change 
which is -0.7608 ± j0.3934 in table 4. The eigenvalue 
associated with the roll mode (0.8625 in table 1) is 
reduced slightly to 0.7975 in table 4. 

The next F&M example illustrates results for a 
double control-effector failure case. At 0.25 second the 
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0 5 10 15 20 


Time, sec. 

Figure 5. Example of system response to an 
aerodynamic change due to missing LE and LSSD 
surfaces using a force and moment approach. 

yaw nozzle goes to a stuck maximum position of 15° 
left, and at one second the left AMT goes to a stuck 
maximum position of 60°. Results from two simulation 
cases are included in figure 6. One example is for the 
case where the control system instantly knows about the 
failures and reconfigures (solid lines) while the second 
example is for the case where the control system does 
not know about the failures (dotted lines). The 



0 5 10 15 20 


Time, sec. 

Figure 6. Double failure case example using a force and 
moment approach. 

command is illustrated by long dashes but the reference 
signal is not included in this figure to avoid clutter. The 
reconfigured example is clearly superior although both 
cases successfully complete the simulation run. The 
non-reconfigured example does a decent tracking job, 
illustrating robustness to these failures. In the 
reconfigured example the right AMT is used to help 
counteract moments produced by the failed left AMT. 

In both cases, the remaining control effectors (not 
shown) go to new positions to compensate for the 
forces and moments generated by the stuck AMT and 
the missing thrust moment from the yaw nozzle. 
However, in the non-reconfigurable example the left 
tail is saturated for a short period during the simulation 
causing the errors seen in the figure. 
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Moment (M) Approach 

As discussed in the Flying Qualities section, the 
main change for the moment approach is that control 
variables are the three angular rates ( Pb and 

commands to the DI loop are their respective angular 
accelerations ( p b des , q bides , and r bdes ). First-order 

command models are used for the rates with 
bandwidths (O iat , co lon , and co dir respectively. An 

outer loop is wrapped around each of the inner loops to 
control slower variables a , (3 , and <p with frequencies 

co a , (Op, and (o phi . 

Table 5 contains the parameter values in the 
following simulations. Variables co equivaient and 

C equivalent are defined in equations 34 and 35. 


Table 5. Command Model Parameters, M Approach 


Parameter 

Value 

^equivalent 

C equivalent 

"to 

2 

1.22 

0.82 

®phi 

0.75 

a lon 

4.48 

2.8 

0.8 

(0 a 

1.75 

®dir 

3.6 

1.8 

1.0 

(Op 

0.9 


Figure 7 shows comparison between a baseline 
simulation and a simulation with three known failures. 
Since slower variables are being controlled in this 
example, longer step responses are used as opposed to 
doublets. In these simulation runs, is commanded to 

5° during the first 5 seconds and then back to null. At 
7.5 seconds into the simulation, a 50° bank-angle step is 
commanded, and between 1 1 and 16 seconds a 5° 

(X -step is commanded. 

The top three plots show the commands (long 
dashed line) and responses (solid line for baseline and 
dashed line for failure case). The remaining 9 plots 
represent control signals as measured at the actuator 
output. Only the leading edge outboard flaps are not 
shown since they are insignificant. Both baseline and 
failure case results are included for each plot. 

Three different failures are demonstrated: 1) at 

0.25 second the yaw nozzle fails at 0°, 2) at 1 second 
the left elevon is missing, and 3) at 1 second the left 
AMT goes to maximum deflection of 60°. The missing 
elevon and stuck AMT cause aerodynamic changes. 

One main conclusion from these results is that the 
reconfigured control system behaves almost as good as 
the original control system. However, the plots 


illustrate that large changes in actuator control positions 
are required to accommodate the three failures. 

Robustness Evaluation 

Stability and robustness evaluation was done 
experimentally using a simulation tool named 
RASCLE 23 (Robustness Analysis for Simulation-based 
Control Law Evaluation) that was recently developed 
under contract to NASA Langley Research Center. 
Basically, RASCLE inserts selected uncertainties 
(gains, delays, bias, noise) into the nonlinear simulation 
and performs a multitude of simulation runs to 
determine the uncertainty values that cause the control 
system to become unstable. In particular, RASCLE 
finds the stability boundary point (or surface for several 
simultaneous uncertainties) within a predefined 
tolerance. For each uncertainty combination the 
nonlinear simulation is first excited by input commands 
(doublets in this paper), and then selected signals are 
monitored for stability. RASCLE performs a power 
spectral density analysis on the monitored data after 
each simulation to determine stability and the dominant 
frequency for each monitored signal. One major 
advantage of using this type of analysis is that 
nonlinear-variables and functions are included. 

RASCLE assumes a stable system with a nominal 
gain uncertainty equal to one and a time delay of zero. 
Gain uncertainties are separated into two parts, gains 
greater than nominal and gains less than nominal. Time 
delays can also be converted to phase angle by using 
the dominant frequency that has been determined. 

For this paper, gain and delay uncertainties were 
inserted separately at several locations to find the 
boundary points. Referencing notation in figure 1, 
uncertainties were inserted at: 1) y and other state 

measurements, 2) accelerometer measurements for x 0 , 
3) actuator control commands 8 cmd , 4) actuator 
position measurements S G , and 5) control derivatives 
B 0 where all derivatives for any control are treated 
together. Specifically, the following measurements 
have been made for both the F&M and the M 
approaches: 

1. y and other state measurements: body rates 

( Pb'Qb' r b )’ linear velocities (u, v, w), angles 

2. Accelerometer measurements for x a (see 
definitions in the Accelerometer Measurements 
section): a x { , a y% [t a n X , a x2 . a yf2 , a n2> a x> 3 , 

a y , 3 » a n . 3 
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SPF, deg. RE, deg. LE, deg. P.deg. < * >,de8 ' a, deg. 



Time, sec. Time, sec. 


Figure 7. Moment-control showing baseline and 3-failure example. 

3. Actuator Control commands S cmd (see definitions 5. Control derivatives B 0 : using columns for 

in the Simulation Results section): LE, RE, SPF, derivatives associated with actuators LE, RE, SPF, 

LAMT, RAMT, LSSD, RSSD, LOLEF, ROLEF, AMT, SSD, LOLEF, ROLEF, PV, and YV. 

PV, YV 

4. Actuator position measurements S Q (same 1 1 
actuator definitions as above but different points) 
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Force and Moment (F&M) Aj 

Dproach 

Moment (M 

) Approacl 

1 

Variable 

Gain > 

Gain < 




Gain < 


Phase 


db 

db 

sec 



db 

sec i 

deg 

state variables r 


Pb 

>20 

<-20 

>2 

N.A. 

>20 

<-20 

0.6094 

86 

% 

12.7 

<-20 

0.1328 

58 

17 

<-20 

0.2969 

89 i 

r b 

>20 

<-20 

0.2656 

121 

>20 

<-20 

0.625 

88 

a 









p 

>20 

<-20 

0.3281 

58 

>20 

<-20 

1.5 

79 

accelerometers , 


a n, 1 

>20 

<-20 n 

0.3906 

76 

>20 

-5.2 

0.2344 

29 

a n, 2 

8.1 

<-20 

0.2281 


2.6 

<-20 

0.1094 

78 

a n,3 

7.4 

<-20 

0.0797 

73 

3.5 

<-20 

0.1562 

85 

actuator control commands — , 


LE 

10.2 

<-20 

0.4219 

163 

12 

<-20 

0.6875 

169 

RE 

12 

<-20 

0.4219 

163 

IEzEHI 



N.A. 

LAMT 

19.1 

<-20 

>2 

N.A. 





RAMT 

>20 

<-20 

>2 

N.A. 


IBISsM 

<-20 

L>1 

control derivative columns 


YV 14.7 

<-20 

>2 

N.A. 

>20 1 < -20 | > 2 

N.A. 


Gain uncertainty boundary limits were set at a 
factor of 10 (± 20 db) and the time delay limit was set 
at 2 seconds. Uncertainty results are shown in table 6. 
Only those variables with at least one uncertainty 
within the boundary limit for which the system is 
unstable are shown in the table. The rest of the 
variables have large uncertainty limits. Phase angle has 
been calculated where applicable; time delays greater 
than the bound do not have a boundary frequency and 
are not applicable (N.A.). Tolerance for the phase 
angle calculations is ±9° based upon the number of 
samples taken for each simulation. All data is for an 
initial trim condition of Mach 0.7 at 25000 feet. It 
appears that the only variables that do not meet military 
specification (6 db and 45°) are the normal 
accelerometers using the M approach. All of the 
normal accelerometers have at least one low gain 
margin and a n { has a low phase margin. A probable 

reason why the F&M approach appears more robust for 
the normal accelerometer uncertainties is that the 
longitudinal DI controller uses a pitch rate command 
whereas in the M approach an a -command is used in 
an outer loop to drive the DI pitch rate command. 
Normal acceleration is used in the nonlinear gain, as 
described in the Flying Qualities section. These results 
indicate a need to integrate a good accelerometer FDI 
(failure detection and isolation) system. 


RASCLE initially found poor margins in the 
actuator control commands for TE and RE. The main 
reason is that the actuator loop is effectively an 
integrator created by feeding back the actuator 
position S 0 (see figure 1). Uncertainties in this loop are 

accumulated and could drive a control into saturation. 

This problem was reduced by adding a simple first- 
order actuator reference model driven by 5 cm ^ . 

Outputs from the reference model are compared with 
corresponding outputs from the real actuator and the 
signal with the lowest absolute value is selected. When 
the real actuator outputs are used, the reference model 
states are updated to the actual states. A reference 
model is used for all actuator signals. 

Conclusion 

A modified version of nonlinear dynamic inversion 
provides the underlying framework of the proposed 
online design control law. The proposed control is less 
dependent on the onboard model than other dynamic 
inversion based methods in that it only requires 
information on control derivatives. No stability 
derivatives are required in the implementation. In 
failure, aerodynamic changes cause force and moment 
changes on the aircraft. The force and moment changes 
can be measured by the accelerometers. The proposed 
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control exploits this result to regain control of the 
damaged vehicle. Sensitivity to the onboard model is 
reduced in favor of obtaining a more robust measure of 
the vehicle’s accelerations. 

Implementation issues associated with obtaining 
the vehicle’s linear and angular accelerations have been 
addressed. Angular acceleration data is derived from a 
unique implementation of linear accelerometer 
measurements, in addition to differentiated rate signals, 
using a least squares solution. Based on the angular 
acceleration, the measured linear accelerations are 
translated back to the eg and averaged. The 
accelerometer implementation used could allow a 
significant amount of sensor redundancy and noise 
reduction, although these aspects have not been 
explored. 

In an attempt to further promote robustness, two 
control structures involving two sets of control 
variables were investigated. The first, the force and 
moment (F&M) approach requires both linear and 
angular control derivatives. The second, moment (M) 
approach requires only the angular control derivatives. 
The derivatives required are determined by the control 
variables specified in the inner dynamic inversion loop. 
The F&M control uses a mixture of force type control 
variables (a and P ) and moment type control 

variables {Pb^Qb< r b )• The M control uses only 
moment type variables in the inner dynamic inversion 
loop. Command models containing flying qualities 
were developed for both approaches, with M control 
requiring a more elaborate outer-loop structure 
involving force type control variables to obtain 
handling comparable to the F&M approach. 

Simulation results indicate both methods are 
capable of online control design subject to multiple 
failures during complex maneuvering. The baseline 
simulation results illustrate excellent tracking for the 
F&M and M control approaches. The F&M 
performance is slightly better than the M performance. 
However, the F&M approach utilizes more involved 
control variables and more control derivatives in its 
inner loop. The M approach utilizes non-linear 
expressions to relate the outer-loop force variables to 
the DI moment variables. Both controls have been 
shown to accommodate aerodynamic changes without 
the use of stability derivative information. The system 
response looks good when there is enough control 
power to compensate for the failures. An F&M 
example illustrated robustness to some failures when 
the failure information was unknown. 

An experimental simulation tool (RASCLE) was 
used to evaluate the control system’s stability 
robustness. Gain and phase uncertainties were inserted 


at five major control system locations and 48 variables 
were evaluated. Only the normal accelerometer 
uncertainties resulted in the control not meeting the 
desired gain margin (6 dB) and phase margin ( 45° ) for 
the uncertainty used. This analysis indicates the need 
for a good FDI system to eliminate defective 
accelerations. 
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